Spin-orbit Coupling Effects on the Superfluidity of Fermi Gas in an Optical Lattice 
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We investigate the superfluidity of attractive Fermi gas in a square optical lattice with spin-orbit 
coupling (SOC). We show that the system displays a variety of new fllling-dependent features. At 
half filling, a quantum phase transition from a semimetal to a superfluid is found for large SOC. 
Close to half filling where the emerging Dirac cones governs the behaviors of the system, SOC tends 
to suppress the BCS superfluidity. Conversely, SOC can significantly enhance both the pairing gap 
and condensate fraction and lead to a new BCS-BEC crossover for small fillings. Moreover, we 
demonstrate that the superfiuid fraction also exhibits many interesting phenomena compared with 
the spin-orbit coupled Fermi gas without lattice. 
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The spin-orbit coupling (SOC) plays a central role in 
the investigation of novel topological states in solid state 
physics [1, 2} • This has stimulated tremendous interests 
in creating artificial non-Abelian gauge fields in ultracold 
atom systems [sj. The successful realization of SOC in 
both Bose-Einstein condensate (BEC) [j, |5[ and Fermi 
gas p, lZI opens up a new avenue towards studying the 
rich physics of spin-orbit (SO) coupled ultracold atoms 
[8|4l3|. One of the important advances is that SOC was 
shown to have fundamental effects on the superfluidity 
of 3D [ii-Eil and 2D [11, |20| continuous Fermi gases. 

On the other hand, the attractive Fermi gas subjected 
to an optical lattice [2l|, |22| has made it possible to simu- 
late the negative- [/ Hubbard model, a basic model for the 
superconductivity of many solid state materials [23| . In 
particular, the on-site attractions can induce deep bound 
states, which cause the conventional BCS-BEC crossover. 
Recently, SOC has been combined to optical lattices for 
repulsive ultracold gases and predicted to lead to many 
interesting phenomena 2J-|26|. Nevertheless, the super- 
fluidity of SO coupled attractive Fermi gas in an optical 
lattice remains a new frontier to be explored. 

In this Letter, we study the Fermi gas subjected to a 
square optical lattice with SOC. Such a system can be de- 
scribed by a generalized negative- /7 Hubbard model. We 
show that, the combination of SOC and lattice can give 
rise to various new features that depend on the fillings. 
Remarkably, there develops a quantum phase transition 
(QPT) from a semimetal to a superfluid for large SOC at 
half filling, with the critical interaction Uc/t ::^ 3.11 (t is 
the hopping amplitude). For close to half filling, we show 
that the emerging Dirac cones governs the behaviors of 
the system, which tends to suppress the BCS superfluid- 
ity. By contrary, SOC can significantly enhance both the 
pairing gap and condensate fraction and lead to a new 
BCS-BEC crossover for small fillings. Compared with 
the SO coupled Fermi gas without lattice, such opposite 
filling-dependent behavior of SOC is rather unique as it 
can only be induced in the lattice system. Furthermore, 



we investigate the superfiuid fraction, which also exhibits 
many unusual characteristics in contrast to the continu- 
ous Fermi gas. 

We consider a system of two-component Fermi gas 
moving in an optical square lattice. In the tight bind- 
ing approximation, the Hamiltonian reads 

<ij> (J(t' 

-U^Ui^nn - fi^rii, (1) 

where t is the overall hopping amplitude and c]^ is the 
creation operator for spin-up (down) fermion a =t,i 
at site i. The nearest sites tunneling matrices Rij = 
g*^(^j-^i) with A = X{axjCry) the non-Abeliangauge- 



field [2J-I2i 

Fig. IHa)). 



A is the strength of Rashba SOC 



(see 



Here, the diagonal term of Rij denotes the 



spin-conserved hopping, while the non-diagonal term can 
be realized by the Raman laser assisted spin-flipped tun- 
neling [29| . U is the on-site attraction strength which can 
be tuned by Feshbach resonances and /i is the chemical 
potential, n = (n^^ + n^^) is the filling factor. 

Fig. mb) shows the band structure of non-interacting 
fermions, where SOC lifts the spin degeneracy and gives 
rise to two split Rashba bands. Remarkably, the two 
bands intersect linearly at F = (0,0), M = (tt, 0), (0,7r) 
and K = (7r,7r). The zero energy Fermi surfaces at half 
filling is shown in Fig. [D^c), where we have a particle 
(hole) Fermi-pocket around F (K) which is associated 
to the up (down) Rashba band respectively. Note that, 
there always exist two zero energy Fermi points at M 
for any A 7^ 0. Specifically, when A = 7r/2, both the 
particle (hole) Fermi-pockets shrink to Fermi points at 
zero energy, and there develops a semimetal with four 
Dirac cones at F, K and M. Fig. [IJd) shows the density 
of states (DOS) p{E) of single-particle excitation over the 
regime A G [0,7r/2], we see that when A = 7r/2, p{E) oc 
|£^| which vanishes linearly around zero energy. 
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FIG. 1: (color online) (a) SO coupled square optical lattice, 
where Rx{y) denotes the non-Abelian hopping matrice along 
x{y) direction, (b) Energy bands of non-interacting fermions 
with A = Stt/IO for illustration, (c) Zero energy Fermi sur- 
faces at half filling, where the filled particle (blue) and hole 
(red) Fermi pockets represent the up/down Rashba band re- 
spectively. The green dots denote two zero energy Fermi 
points at M. (d) Density of states p{E) over the regime 
AG [0,7r/2]. 



We start by writing the partition function in the 
imaginary-time path integral Z = JP[?/5, ?/^]e~'^[^'^], 

where S[i),^lj] = Jq dr[}2^i^adri^a + i^(V^,V^)] with 
tjj = (t/j^^ipi)^ representing the Grassmann field vari- 
ables. Then, by decoupling the attractive term of 
Eq. ([1]) in normal and anomalous channels through 
a pairing field Ai{r) = UilJii{r)ilJi^{r) and introduc- 
ing (V^kt^V^ki, V^-kt? V^-ki)"^, we obtain the effective ac- 
tion after integrating out the fermionic field 6'eff = 
E^ lo dr^-^^ - ^Tr In G'^ + /3 Ek ^k- Here we have 
ignored the constant term NUn? /A {N is the number of 
lattice sites) and the inverse Green function is given by 



dr + e_k + ^k 

iayAi{r) 



-iayAi{r)_ 

9r - £k + Ak 



(2) 



with ^k = — 2tcosA(cos/caj + cos ky) — fl and Ak = 
— 2tsinA(sinA:a^cra^ + sin/c^cr^), where ft = /j. -\- Un/2 
is the scaled chemical potential. Furthermore, we set 
Ai(r) = A + (5A and write G'^ = G"^ + S with 
G~^ = 5~^|Ai(r)=A- Then, the effective action can be 
expanded to the second order of fluctuation S as ^eff — 
So^AS with So = ^E |Ap + i 5:^^^^^[f (sk-^k,.)- 
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)], and AS = ^J-Hq)SA{-q)5A{q) 



jjj:^5A(-q)5A{q) + iTr[G(fc)S(-9)G(fc 
Here k = (k, zk;^), q = (q, iz^n), and £^k,± = 



iMq)]- 
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At the mean-field level, the 



many-body ground state of the system can be derived 
by minimizing So/{Nl3) with respect to A and /i, and we 
have the following gap and Fermi density equations 



k,i/=± ' 



pEk,. 



), 



- y 



^k , 1 //^£^k,ryx 

,. / tannf —). 

k,z/=± ' 
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Before proceeding, it's useful to consider the large at- 
traction limit with U jt ^1. In this case, the stan- 
dard degenerate perturbation theory can be applied for 
Eq. ([1]) through the canonical transformation q^ -^ q^ 
and Ci^ -^ (-l)*^+*^cj^ [3^. For any band filling, we 



<i3> 



with ^k,± = £k±2t sinA/C being the two Rashba branches. 



can derive an effective spin model i^spin — J^ 
^3 ~ 2/i ^^ SI in the presence of arbitrary SOC. Here 
J = At^ /U and the pairing field operator becomes the 
transverse magnetic operator. For /i 7^ 0, the antiferro- 
magnetic order in XY plane is equivalent to the pairing 
order of Eq. ([1]). Therefore, we conclude that SOC does 
not have any influence on the superfluidity of the system 
in the large U limit, where all the fermionic atoms form 
tightly bound molecules and give rise to a Kosterlitz- 
Thouless transition of BEG [31[. In the following, we 
shall focus on the more interested weak and intermediate 
attraction regions. 

The pairing gaps at zero temperature are illustrated in 
Fig. El First for half filling (n = 1), Fig. [2](a) shows that 
the BGS gap decreases monotonically with respect to A. 
This could be understood that, the Fermi pockets around 
F and K (see Fig. HJc)) tend to form the Fermi points 
at Ef = by increasing SOG, which causes a suppres- 
sion of DOS at zero energy (see Fig. [Hd)). Specifically, 
when A = 7r/2 the system becomes a semimetal, which 
is expected to be stable towards small attractions. On 
the other hand, when U /t ^ 1 the system should sup- 
port a superfluid state of bound molecules as indicated 
by the effective spin model i^spin [32|. Hence, there must 
undergo a significant QPT from a semimetal to a super- 
fluid by increasing attractions, see the thick vertical line 
of Fig. [2] (a). In the inset, we show that the critical value 
[/c/^ — 3.11, above which a finite gap develops. 

Such scenario would be affected remarkably by dop- 
ings as shown in Fig. [2] (b). Without loss of generality. 
We focus on the hole doping case due to the particle-hole 
symmetry of the system. First for small dopings, simi- 
lar to that of half filling, the superfluidity is governed by 
emerging Dirac cones at zero energy and A is suppressed 
by increasing A, see n = 0.95. However, when close to 
A = 7r/2, the doping would make the QPT at half filling 
unstable and opens a gap. This produces a nonmonotonic 
behavior of A with a minimum at Amin- While for large 
dopings, the situation is entirely changed, where the in- 
fluence of Dirac cones would diminish and SOG induces 
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FIG. 2: (color online) (a) Plot of A versus A at half filling for 
different IJ ji. Inset shows a semimetal-superfluid QPT with 
[/c/t 2^ 3.11. (b) Plot of A and (c) [i versus A with V ji = 2 
for different fillings, (d) Pairing gap as a function of JJ ji for 
two typical fillings with n = 0.1 (red) and 0.95 (blue). Solid 
and dashed lines represent A = and 7r/2 respectively. 



a new BCS-BEC crossover with A being significantly en- 
hanced, see the dash dotted line of n = 0.1. This can be 
understood by solving the two-body problem of Eq. ([T]), 
which is determined by V~^(iVn ^ a; + z0+,q = 0) = 
as cj + 2/i = —Eb and we arrive at 



TT ON /-^ 



1 



u 



2^ k^± 2(a, 



Eq) — Eb 



(4) 



where Eq denotes the lowest energy of ^k,- branch. In the 
absence of SOC (A = 0), the binding energy \EB\/t ~ 
in the weak attraction region U/zt < 1 (z = 4 is the 
number of the nearest neighbor) and becomes very large 
for U/zt :^ 1, evolving from loosely local pairs (BCS) 
to tightly bound molecules (BEG) [23|. However, when 
SOG is added to the lattice, \EB\/t will be significantly 
enhanced (left panel of Fig. [3]) due to the increasing of 
DOS around Eq and effective shrinking of the bandwidth 
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FIG. 3: (color online) Binding energy EB/t as a function of 
U/t (left panel) and SOG strength A (right panel). 

(see Fig. md)) [33|. In particular, the right panel shows 



a remarkable grow of \EB\/t from nearly zero in the weak 
attraction regions, which signifies the formation of SOG 
induced bound states, see U/t = 2,3 for example. 

In general, such filling-dependent effects arise from the 
unique features of combination of SOG and lattice. The 
opposite behaviors of A versus SOG upon dopings indi- 
cate that the system evolves from the Dirac cone domi- 
nated physics near half filling to the SOG induced BGS- 
BEG crossover at small fillings. In this respect, the gap 
behavior at small fillings is reminiscent of SO enhanced 
pairing in the unitary Fermi gas 14I42QJ. However, we 
note that the lattice also plays nontrivial roles even in 
this limit. In Figl2](c), we see that /i increases with SOG 
at n = 0.1, which differs from the familiar results in the 
continuous system where /i is decreased by SOG. 

The opposite roles of SOG in lattice can be clearly seen 
in Fig. [2] (d), where we plot A as a function of U/t for 
two typical fillings. We take A = 7r/2 for illustration 
and show that the strong SOG can remarkably enhance 
(n = 0.1) or suppress (n = 0.95) the pairing gaps of the 
conventional BGS-BEG crossover of negative- [/ Hubbard 
model, especially in the weak and intermediate attraction 
regions. While in the large attraction limit, A versus 
U/t approaches the A = results, which means the SOG 
effects diminish according to the effective theory of i^gpin- 

To get more insight of the unusual properties of 
this system, we now explore the condensate and 
superfluid fractions. First, the condensate density 
^c = ^Ek,a,a' KV^kaV^-ka')P H , whcrc the sin- 
glet and induced triplet pairing fields (V^ktV^-kt) = 



.A„-i^k " 



^ E. ^ and (V^ktV^-ki) = -f E. V^k,. with 
^k = SiTg{smkx + i sin ky). While for the superfluid 
density, we impose a phase twist on order parameter 
A -^ Ae*^^"^^ by a local unitary transformation ipj -^ 
^^.g*6'(rj)^ Then, the inverse Green function can be writ- 
ten as G-i[A, V(9] = G-i[A] + S[V(9]. After lengthy but 
straightforward calculations, we derive a classical phase 
variation model H = \j j d^Y[{d:^Of^{dyef] with J the 
phase stiffness. Therefore, the superfluid density can be 
defined as ps = oiiv' which reads 



cosAv^ , sinAv-^ i/^k,i/ svo? kx ^ , ,pEvy^ 

^ -tanh( — —^) 



-2^cosA;^nk + -r^2^- 



N ^ - -' N ^ 2^k,. /C " ^ 2 



2t_ 
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y^ f'{E\<i,v) sin^ kx f cos A - 



sin A cos kx 



sin A ^^ e\ + v2t sinA /C^k + A^ sin ky cos^ k^ 



E' 



k,i/ ' 

X tanh( —). 



K? 



(5) 



r^E^^ 



Here nk = 1 — Ezy=± 2e' tQ-^h( ^'"^ ) and the third 
term vanishes at T = 0. Not that, although the first 
and fourth terms bear similarities with the continuous 




FIG. 4: (color online) (a) Condensate density ric and (b) su- 
perfluid density ps (both divided by n) at T = as a function 
of SOC for different fillings, where we take U/t = 2. (c) Plot 
of Tic and (d) ps versus U/t for two typical fillings with n = 0.1 
(red) and 0.95 (blue). Solid and dashed lines represent A = 
and 7r/2 respectively. 



dependent features. At half filling, we find a QPT from a 
semimetal to a superfluid for large SOC. While upon dop- 
ings, the system evolves from the Dirac cone dominated 
physics near half filling to the SOC induced BCS-BEC 
crossover at small fillings. Moreover, we show that all 
the pairing gap, condensate and superfluid fractions ex- 
hibit many interesting physics, which differ qualitatively 
from the SO coupled Fermi gas without lattice and the 
conventional negative- /7 Hubbard model without SOC. 
We hope that this work will trigger new exciting inter- 
ests to the SO coupled optical lattice physics, and may be 
useful for the study of superconductivity of future solid 
state materials with SOC. 
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system jl9|, |20|, they may become negative in the lattice. 
Fortunately, the new second term which exists only in the 
SO coupled lattice system, can stabilize the superfluidity 
and gives rise to many intriguing features. 

Fig. |4] (a) shows the condensate fraction ric versus 
SOC with the evolution of fillings. For n = 0.1, we see 
that Tic increase dramatically above a characteristic value 
Ac, which agrees with the formation of the SOC induced 
bound states. The characteristic value Ac grows with 
increasing fillings, and until n o:^ 0.7, ric begins to de- 
crease with respect to A. This indicates that the BCS 
superfluidity would be suppressed at large fillings, see 
n = 0.95 for example. On the contrary, the superfluid 
fraction ps is always suppressed by SOC and generally 
decreased with increasing n, as shown in Fig. [41(b). Sig- 
nificantly, there exhibits a characteristic minimum of A, 
which moves rightwards when n is increased. 

In Fig. [4] (c) and (d), we plot Uc and ps with respect 
to U/t for two typical fillings. Qualitatively different 
from the conventional BCS-BEC crossover, the conden- 
sate fraction is remarkably enhanced (n = 0.1) or sup- 
pressed (n = 0.95) by the strong SOC. Conversely, the 
superfluid fraction is destroyed by both increasing U/t 
and A. Inter set ingly, ps exhibits quite a nontrivial behav- 
ior at large fillings. The presence of strong SOC can dra- 
matically suppress ps for the weak attractions and causes 
a maximum of ps located at the intermediate crossover 
region, see n = 0.95 for example. While in the large at- 
traction limit, both Uc and ps will approach the results 
without SOC. 

In summary, we have shown that the SO coupled Fermi 
gas in an optical lattice displays various new filling- 
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